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INTRODUCTION 


This  is  the  final  report  for  Contract  N00014-81- 
C-2355.  It  covers  approximately  40  man-days  of  research  per- 

i 

formed  by  ^erkeley  Scholars,  Inc.  from  April,  1981  through 
April,  1982. The  research  involved  the  theoretical  analysis 


of  charged  particle  beam  and  other  plasma  systems. 

The  major  emphasis  of  the  research  was  on  modeling 
the  free-electron- laser  (FEL)  oscillator.  The  approach  was  to 
utilize  a  first  order  description  for  the  evolution  of  rela¬ 
tivistic  electron  beam  pulses  in  a  cavity  with  an  applied 
wiggler  magnetic  field  to  develop  a  model  for  the  modal  elec¬ 
tromagnetic  energy  in  a  FEL  oscillator.  In  the  following 
TECHNICAL  DISCUSSION  we  present  the  de^Lls  of  our  research. 


technical  discussion 


During  the  performance  period,  Berkeley  Scholars, 
Inc.  was  funded  for  approximately  40  man-days  to  examine 
several  areas  of  charged-particle  beam  research  and  other 
plasma  systems.  These  included  relativistic  charged-particle 
beam  propagation  and  the  free-electron-laser .  The  free-elec- 
tron-laser  is  of  great  interest  due  to  its  potentially  high 
efficiency  (-20%) ,  high  power  operation  (few  MW) ,  and  continu¬ 
ous  frequency  tunability.  The  major  emphasis  of  our  research 
was  on  modeling  the  free-electron-laser  (FEL)  oscillator. 
Discussion  and  details  of  our  research  follows: 

The  FEL  Oscillator 

Figure  1  is  a  pictorial  description  of  a  free- 

electron-laser  at  t  =  0.  Nq  is  the  average  density  of  the 

incoming  electron  beam  pulses  of  width  L^  and  separation  Lc. 

The  cavity  mirrors  are  at  z  =  0  and  z  =  L  with  the  wiggler 

magnetic  field  applied  between  z  =  L^  and  z  =  L  +  L  .  For 

r  o  o  w 

a  relativistic  electron  the  equation  of  motion  is 

St  {m^  =  _e  (k  + 


where 


Introducing  the  vector  potential  A  and  the  scalar 
potential  $  such  that 

B  =  V  x  a  and  E  =  -V$-^ 

—  —  —  c  at 

and  assuming  for  the  canonical  model  of  the  FEL 


A  =  A  (z, t) e  +  A  (z,t)e 
—  a  *  y  -y 

4>  =  $(z,t) 

yields 

mYVj.  =  ^  where  Vj.  »  vxex  +  vyey 
and  in  the  tenuous  beam  limit  where  *  is  negligible 

A2 

T  * 


+  z 

1 

bo  =  -JP  (®xcoskwz  +  eysinkwz)  H(z-L0)H(L0+lw-z) 

where  Bw  and  kw  are  the  amplitude  and  wavevector  of  the  wiggler 
field  and  H  is  a  heaviside  function  yields 

v.  "  "(s#£ )( ^  &  +  •  (i 


z  ImcY y  l  c2  3t  3z/ 


Linearizing  to  first  order  with 


A  -  A  +  A 
—  -o  —i 


z  =  z_  +  v  t 
o  o 
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cavity 


Expanding 


in  terms  of  the  vacuum  modes  of  the 


A 
— 1 


-iw  t 

a  (t)e  sink  z 

— n  n 


where  a  =  a*,  w  =  and  k 

— n  — n  n  L  n 


n  it 


L 


(2) 


and  substituting  equation  (2)  into  equation  (1)  yields  for 

L  <z<L  +  L 
o  o  w 


•a  e 
— n 


i(k„+kw> 


zo+il,nt 


'^ne 


-i(kn+kw>  z'>-ivnt 


where  yn  =  <kn+kw) v0~  I wnl  and  lin^n^*1*  since  Z1  and  z, 
vanish  if  the  electron  in  question  has  not  entered  the  wiggler 
z^t)  =  /  dxtt-Tj’z^  (t)H(z0+v0t-L0)H(L0+Lw-Z0-v0t) 

and  it  follows  that 

L  -z 

for  t  <  z  (t)  =  0  , 

vo  1 

L  -Z_  L  +L  -Z 


The  above  calculations  are  for  an  individual  electron 
within  the  beam  pulse.  Therefore  the  beam  density  is  defined  by 

N  =  i  ^  6  (z-z j  (t) ) 

j 

where  o  denotes  the  cross-sectional  area  of  the  electron  beam. 
The  average  density  in  the  absence  of  the  high  frequency  field 
is  then 

N  =  4  <  ^  <5(z  -  z  .  -  v  t)  > 
o  o  ^  03  o 

j 

where  <...>  denotes  the  ensemble  average.  Consequently,  the 
equation  for  the  first  order  vector  potential  is 


i  (k_+k..)  z+iw_t 


and  k>  =  b  for  n>l  . 
— n  — n 


Equations  (6)  and  (7)  suggest  writing 


2m  -  <£x  +  “me 


2n  =  <£x  +  iSy*  6ne 


-iy  t 


Forming  the  scalar  product  of  equation  (6)  with  (ex  - 


yields 


a  =8  +  M  a 

m  Mn  mn  n 


or  in  matrix  notation 


a  =  8  +  Ma 


where 


L  +L  -v  t 
*  o  w  o 


M  +  iu  6 
mn  m  mn 


/  eB  \2  2  f  °  w  u 
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The  objects  of  physical  interest  are  the  <[cij|2> 


which  measure  the  energy  per  line. 


Defining  the  matrix 


C  =  <aa> 


where  ~  denotes  the  Hermitian  conjugate  we  can  then  formally 
solve  equation  (9)  for  C. 


a(t)  =  J  dTP(t,T) &(x) 


where  P(t,x)  =  Q  (t)Q(x),  Q  =  -QM,  Q(0)  =  I,  and  I  is  the 
identity  matrix.  Note  that 

P(t,x)  =  M(t) P (t, x)  and  P(t,t)  =  I  . 

Therefore 

C  (t)  *  f  dx  /*.  dx'P(t,x)<e(x)3(T')  P  ( t ,  x  ' ) 

Jo  Jo 


and 


where 


C(t)  =  MC  +  CM  +  D 


dx<3(t) 8(x) >P(t,x)  +  /  dxP (t,x) <8(x) 8(t) >  .  (11) 


Clearly 


TRACE  C  =  J-  (ensemble  average  electromagnetic  energy) 


* 


IH 


.  ^ 


.  j 


and 


(TRACE  C)  =  TRACE  (M+M) C  +  TRACE  D  . 


(12) 


Consequently  the  Hermitian  par-  of  M  determines  the  total  gain. 


From  equations  (7)  and  (8) 


a  a)  to 


rYok„L , 
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)  <»n(t)l5m(T)> 
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(13) 
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For  t>T  equation  (13)  involves  at  most  one  pulse  and  requires 
that  the  pulse  be  in  the  wiggler  in  order  to  emit.  From 
equations  (11)  and  (13)  when  (k^-k^) L.<1 


&  -i(Vk„><iLc4Lb> 


NL,  e 
b 


H{L  +L  + jL  -V  t) 
o  w  J  c  o 


(14) 


where 


NL, 


r~iL c 

=  /  dz  N  <z. 


0) 


and  D  is  a  positive  definite  matrix.  Therefore  from  equations 
(10) ,  (12) ,  and  (14)  the  evolution  of  the  ensemble  average  of 
the  electromagnetic  energy  may  be  calculated. 
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Figure  1.  Pictorial  description  of  a  free-electron-laser  at 


